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GeV gamma-rays and TeV neutrinos from very massive 
compact binary systems: The case of WR 20a 
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ABSTRACT 

Massive Wolf-Rayet stars in a compact binary system are characterized by very strong 
winds which collide creating a shock wave. If the wind nuclei, from helium up to oxy- 
gen, accelerated at the shock can reach large enough energies, they suffer disintegration 
in collisions with soft thermal radiation from the massive stars injecting relativistic 
protons and neutrons. Protons collide with the matter of the wind and a fraction of 
neutrons collide with the massive stars producing 7-rays and neutrinos via decay of 
pions. We calculate 7-ray fluxes from inverse Compton e pair cascades initiated by 
primary 7-rays and leptons, which are produced by protons, and the neutrino fluxes, 
produced by protons and neutrons, for the example compact massive binary WR 
20a. From normalization of these cascade 7-ray spectra to the fluxes of the EGRET 
sources, 2EG J1021-5835 and 2EG J1049-5847 observed in the direction of WR 20a, 
we conclude that this massive binary can be detected by the IceCube type neutrino 
detector. The estimated muon neutrino event rate inside the 1 km 2 detector during 1 
yr is between a few up to a few tens, depending on the spectrum of accelerated nuclei. 

Key words: binaries: close - stars: WR 20a - radiation mechanisms: non-thermal - 
gamma-rays: neutrinos: theory - sources: 2EG J1021-5835; 2EG J1049-5847 



The observation of neutrinos from discrete sources 
should give a unique information on the acceleration of 
hadrons in our Galaxy and on the origin of cosmic rays 
at energies above ~ 10 12 eV. Possible production of neu- 
trinos by relativistic hadrons accelerated inside the binary 
systems has been discussed in the case of accreting neutron 
stars (Gaisser & Stanev 1985, Kolb, Turner & Walker 1985, 
Berezinsky, Castagnoli & Galeotti 1985, Anchordoqui et 
al. 2003), energetic pulsars (Harding & Gaisser 1990), 
and accreting solar mass black holes called microquasars 
(Levinson & Waxman 2001, Distefano et al. 2002, Romero 
et al. 2003). For recent review of these models see e.g. 
Bednarek, Burgio & Montaruli (2005) . Most of these models 
predict neutrino fluxes which should be easily detectable by 
the 1 km 2 neutrino detector of the IceCube type and in some 
cases by detectors of the present AMANDA II type. Recent 
searches of the discrete neutrino sources by the AMANDA II 
detector, which base on the analysis of the three-year exper- 
imental data (Ahrens et al. 2004, Ackermann et al. 2005), 
show no clear evidence of a significant neutrino flux excess 
above the background. 

In this paper we show that also compact binary sys- 
tems containing two massive stars can be interesting sources 



of TeV neutrinos likely to be observed by large scale neu- 
trino telescopes. In fact, between 227 Wolf-Rayet (WR) stars 
(van der Hucht 2001), 39% form binary systems (including 
probable binaries), and about 15% of them form compact 
massive binaries with orbital periods < 10 days (see Table 
15 in van der Hucht 2001). One extreme example of a bi- 
nary system, WR 20a, contains two WR stars with masses 
~ 70M Q , surface temperature T = 10 5 T 5 K « 4 x 10 4 K, 
and radii ~ 20Rq (where Rq — 7 x 10 10 cm is the radius of 
the Sun), on an orbit with the semimajor axis a/2 ~ 26i?e 
(Rauw et al. 2004, Bonanos et al. 2004). Let's consider such 
a very massive binary composed of two massive WR type 
stars characterized by very strong winds, with mass loss rate 
M = W~ 5 M- 5 M yr _1 ~ (0.8 - 8) x 1O" 5 M yr -1 , wind 
terminal velocities = 3 x 10 3 W3 km s _1 ~ (1 — 5) x 10 3 km 
s _1 . The radiatively driven stellar winds accelerate from the 
surface of the star according to 113(7") = vjj + (1)3° — v^)(l — 
1/r) 13 w «3°(1 — 1/?"), where V3(r) is the wind velocity at 
the distance, R = tRwr, from the star, v® — O.OIV3 , /3 — 1 
(e.g. Castor & Lamers 1979). The surface magnetic fields 
of WR stars can grow up to B B ~ 10 4 Gs (Maheswaran & 
Cassinelli 1994). For the massive stars in WR 20a binary 
system we apply the typical values of B — 1Q A B3 G, i.e. sig- 
nificantly below 10 4 G. Such stellar winds collide creating 
a double shock structure separated by the contact disconti- 
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nuity which distance from the stars depends on their wind 
parameters. 

It is likely that the binary system WR 20a is responsible 
for one of the EGRET sources observed in this region, i.e. 
2EG J1021-5835 and GeV J1025-5809 or 2EG J1049-5847 
and GeV J1047-5840 (see Thompson et al. 1995, Lamb & 
MaComb 1997). They have flat spectra with the index close 
to 2 (Merck et al. 1996) and fluxes at the level of ~ 10~ 7 
cm -2 s" 1 above 1 GeV (Lamb & MaComb 1997). 

Nuclei present in the stellar winds can be acceler- 
ated at such shock structure either due to the magnetic 
field reconnection (we call this possibility as model I) 
or the diffusive shock acceleration mechanism (model II) 
(see Fig. 0. The acceleration occurs in the fast reconnec- 
tion mode provided that following condition is fulfilled, 



2m e /Am p « W~ 3 /A < 8tt pk B T sh / B 2 h « 3.3 x 



10 Ms/B 3 r V3(r) (although the fast reconnection can 
also occur in specific situations when this condition is not 
met, e.g. Hanasz & Lesch 2003), assuming typical tempera- 
ture of the plasma at the shock T s h = 10 7 K (e.g. Stephens 
et al. 1992), density of the wind p = M/4nR 2 v w w 5.4 x 
10 10 M-^/r 2 v 3 {r) cm -3 . Then, the reconnection occurs with 
the speed close to the Alfven velocity va = B s h/\/4^™pP ~ 
7 x 10 s B 3 vl /2 (r)/rM% 2 cm s . Nuclei can reach maximum 
Lorentz factors of the order of 

^B 2 L 12 v* /2 (r) 



ZeB s hLrecVA 

Am p c 



2.6 x 10° 



r 3 M 



1/2 



(1) 



where A/Z = 2, m p « m n is the nucleon mass, L rcc = 
10 12 Li2 cm is the length of the reconnection region, c is the 
velocity of light, and e is the proton charge. In order to es- 
timate the magnetic field strength after the shock, we apply 
the model for external magnetic field structure in the case 
of strong outflowing gas (Usov & Melrose 1992). According 
to this model the magnetic field is a dipole type below the 
Alfven radius, Ra, becomes radial above R\, and at the 
largest distances, determined by the rotation velocity of the 
star, i.e. at above ~ 10-Rwr, the toroidal component dom- 
inates. We are interested mainly in the region in which the 
dipole and radial components dominate. Above ~ IORwr 
the product of the magnetic field strength at the shock and 
the length of the reconnection region is roughly constant 
and the maximum energy of accelerated nuclei does not de- 
pend on the distance from the star. The Alfven radius for 
the example parameters considered in this paper, vf = 1, 
B 3 = 1, and Ms = 3, is at Ra ~ 1.3-Rwr- Therefore, we 
estimate the magnetic field at the shock in the region of 
radial magnetic field from 

B sh (r) w B s (i? W R/i? A ) 3 (i?A/iish) 2 « 750B 3 /r 2 Gs, (2) 

where R s h is the distance from the center of the star to the 
acceleration region at the shock, i.e. -R s h = rRwK- 

Nuclei can be also accelerated diffusively by the I or- 
der Fermi shock acceleration mechanism. In this case they 
obtain a power law spectrum. Nuclei are accelerated to 7a 
during the time estimated by r aC c ~ (-Rl/c)(c/u w ) 2 , where 
R L = Am p "f A /(eZB sh ) w 8 x 10 9 7 6 r 2 /B 3 cm is the Lar- 
mor radius required to complete acceleration to 7A = 10 6 76. 
The maximum energy of nuclei is not limited by any ra- 
diation loss process but by their escape from the acceler- 
ation region. We identify this escape mechanism with the 
convection of relativistic nuclei with the wind flow along 




Figure 1. Schematic picture of a compact binary of two massive 
stars (WR+WR or WR+O) which create a strong shock as a 
result of their wind collisions. Nuclei are accelerated at the shock 
in the reconnection regions or by the diffusive shock acceleration 
mechanism to the Lorentz factors 7 max - Only nuclei accelerated 
within a part of the shock (marked by thick line within the solid 
angle AC) are accelerated to the Lorentz factors above 7 m i n which 
allow their efficient photo-disintegration in collisions with thermal 
photons from the massive stars. Significant fraction of neutrons 
dissolved from nuclei move toward the massive stars (due to the 
head on collisions) and interact with the matter of the stellar 
atmospheres producing neutrinos. 

the surface of the shock. The escape time is estimated as 
r c « 3R B h/v„ — 10 4 ri?i2/«3(r) s, where the radius of the 
WR star, 7?wr = 10 12 J?i2 cm. The maximum energies of 
nuclei are determined by the condition t c = race Therefore, 
we expect the cut-off in the power law spectrum of nuclei at 
energies corresponding to 

3eZB s hR s hV^ 



7m ax 



cAm„ 



5.3 x 10 J R 12 rv 3 (r)B sh (r). 



(3) 



where B B h is given by Eq. (2). 

The WR stars are the late stage of evolution of super- 
massive stars which have already lost their hydrogen en- 
velopes. Therefore, their winds are mainly composed of nu- 
clei heavier than hydrogen. The heavy nuclei accelerated 
to Lorentz factors described by Eq. Q and Eq. El lose nu- 
cleons due to the photo-disintegration process in collisions 
with the thermal photons from the massive stars if the 
energies of photons in the reference frame of the nuclei, 
E 1 = 3fc B T7A(l + cos6») = 25T 5 7 6 (l + cos<9) MeV, are above 
~ 2 MeV, where fee is the Boltzmann constant, 9 is the an- 
gle between photon and nucleus. The above condition gives 
the lower limit on the Lorentz factor of accelerated nuclei, 



In 



8 x 10 4 /T 5 (1 + . 



(4) 



Most efficient disintegrations occur when photon energies 
correspond to the maximum in the photo-disintegration 
cross-section of nuclei, which is at ~ 20 — 30 MeV with 
the half width of ~ 10 MeV for nuclei with the mass num- 
bers between helium and oxygen (main composition of the 
surface of the WR stars) . Let us estimate the efficiency of 
the photo-disintegration process of nuclei in the considered 
here scenario. The average density of thermal photons from 
both stars (with similar temperatures), at the location of 
the shock, can be approximated by 



liWR 



4as B r 3 ( Rwn Y 
3cfc B V R B h J 



2 x 10 16 -% ph. cm" 3 . 



(5) 



where <tsb is the Stefan-Boltzmann constant. Then, the 
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mean free path for dissociation of a single nucleon from a 
nucleus is Aa 7 = (nwnVA-/) -1 ~ 3.5 x 10 10 r 2 /(AT 5 3 ) cm, 
where for the photo-disintegration cross section the peak in 
the giant resonance is applied, <ta 7 = 1.45 x 10~ 27 A cm 2 
(Karakula & Tkaczyk 1993). For the considered here binary 
system WR 20a, for which R Bh = a/2 and so r = 1.3, the 
characteristic photo-disintegration time scale of nuclei with 
the mass number A, ta 7 = Aa 7 /c « 30/^4 s, is much smaller 
than the acceleration time and the convection escape time 
of nuclei from the acceleration region at the shock (esti- 
mated above). Nuclei with the initial mass numbers between 
4 (helium) and 16 (oxygen) , accelerated to energies given by 
Eqs. and [3] should suffer complete fragmentation. Note 
that nuclei, and protons from their disintegration, can not 
escape far away from the shock region. For typical values 
of the binary system WR 20a, their Larmor radii, i?L, are 
significantly shorter than the characteristic distance scale of 
the considered picture defined by the semimajor axis of the 
binary system. Significant fraction of neutrons from disin- 
tegration of nuclei move toward the surface of the massive 
stars since the probability of dissociation of a single nucleon 
is the highest for the head on collisions of nuclei with ther- 
mal photons. These neutrons propagate along the straight 
lines and interact with the matter of stellar atmospheres. 
On the other hand, protons from disintegration of nuclei are 
convected outside the binary system along the surface of the 
shock structure in a relatively dense stellar winds. 

The relativistic nuclei accelerated at the shock take a 
fraction, £, of the kinetic power of the two stellar winds, 

Pa = tMvl « 6 x 10 37 £M- 5 v 2 (r) erg s" 1 . (6) 

The parameter £ is determined by the efficiency of accel- 
eration of nuclei and by the solid angle AQ. subtended by 
the active part of the shock (see Fig. 1). The active part 
of the shock is determined by the power of the stellar wind 
which falls onto the shock region in which nuclei can be 
accelerated above 7 m i n (see Eq. |1J. AQ. can be estimated 
from the condition of acceleration of nuclei above 7 m i n by 
applying Eqs. (for model I) or [3] (for model II) and Eq.|5| 
By comparing Eq. with Eq. [I] (for 8 — 0°), we estimate 
the maximum distance of the active shock from the star 
on i? m ax ~ 4(_B|T 5 Li2i'3' /2 (r)) 1/ ' 3 _RwR. Similar procedure, 
comparison of Eq. [3]witli Eq. [I] allows us to estimate R^ a3t 
in model II, P max « 10 2 (Ri 2 B 3 T 5 V3(r))R WR . A part of the 
shock solid angle at which nuclei are accelerated above 7111111 
is then AQ = 27r(l — R a / Rmmx) / in , where the distance of 
the shock in the plane of the binary system is approximated 
by R « a/2. For the considered above parameters of WR 
20a, B 3 = 1, T 5 = 0.4, w|° = 1, L 12 = 0.3a/(10 12 cm), 
a/2 = 1.3-Rwr, and Rwr = 20R Q , AQ ~ 0.2 and 0.5 for 
the models I and II, respectively. 

Provided that complete disintegration of nuclei occurs, 
as envisaged above in this scenario, the flux of nucleons dis- 
solved from nuclei in the case of monoenergetic acceleration 
in the reconnection regions (model I) is 

JV„ = PA/m p7A « 4 x 10 34 £M- 5 U3M/76 N s" 1 , (7) 

In the case of nuclei with the power law spectrum (model 
II) the differential spectrum of nucleons is, 

dN n /d-y n dt = Kj~ s , (8) 

in the range 72 = 7 max (Eq. EJ and applying 71 = 10 1 , 



5 is the spectral index, and the normalization constant 
is equal to K = APA/(Arn p In (72/71)) for S — 2 and 
K = AP A (5"2)/(Am p ('yl~ S) -^- s )) for S > 2. The spec- 
tral index of particles accelerated at the shock depends on 
the Alfvenic Mach number of the plasma flowing through 
the shock. It is defined as M — v w /va (see e.g. Schlick- 
eiser 2002), where da is the Alfven velocity defined above. 
For the parameters of WR 20a, we estimate the Alfvenic 
Mach number as a function of distance from the star on 
M ~ r, where the density of the wind and the magnetic 
field strength at the shock location are calculated above. 
We conclude that the outer part of the shock, r> 1, fulfills 
the condition of the strong shock, M 2> 1, although in the 
inner part, the Alfven shock is probably created. Therefore, 
we consider the power law spectrum of accelerated nuclei 
with indexes between 2 (characteristic for a strong shock) 
and 2.3 (for a weaker shock). 

Since the winds of the massive stars are very dense, 
protons from disintegration of nuclei have chance to interact 
with the matter of the winds. The characteristic time scale 
for collision of relativistic protons with the matter of the 
wind is t pp ~ (co-ppp) -1 « 10 4 r 2 V3(r)/M_5 s, where p is the 
density of the stellar wind estimated above. For the applied 
above parameters of the wind, vf = 1 and M-5 = 3, r pp is 
shorter than escape time scale of protons along the shock, 
t c , at distances less than r ~ 5, i.e. in the main part of the 
active shock. 

We calculate the 7-ray spectra from decay of pions 
which are produced by protons in collisions with the matter 
of the wind by integrating the injection spectra of protons 
(monoenergetic or the power law) over the active part of 
the shock and applying the scale break model for hadronic 
interactions developed by Wdowczyk & Wolfendale (1987) 
which is suitable for the considered energies of relativistic 
protons. Only single interaction of proton with the matter 
has been included. These high energy 7-rays originate rela- 
tively close to the surface of the massive stars, r ~ 1.3 — 5, 
and therefore they can suffer absorption in the thermal ra- 
diation coming from the stellar surfaces. In Fig. 2a, we show 
the optical depths for 7-rays as a function of their energies. 
It is clear that most of the 7-rays with energies above a 
few tens GeV are absorbed and initiate the inverse Comp- 
ton (IC) e pair cascades in the radiation of the massive 
stars. In order to determine the final 7-ray spectra which 
escape from the binary system, we apply the Monte Carlo 
code developed by Bednarek (2000) which follow the IC e ± 
pair cascade in the anisotropic radiation field of the massive 
star. It is assumed that primary 7-rays are produced by pro- 
tons isotropically at the active part of the shock. The phase 
averaged 7-ray spectra escaping from the binary system are 
shown in Fig. 2b for the monoenergetic and the power law 
spectra of injected protons, assuming the distance to the WR 
20a equal to 5 kpc (the average value estimated from the 
reported range (2.5-8) kpc, Churchwell et al. 2004). These 
spectra have been normalized to the 7-ray fluxes observed 
from the EGRET sources, 2EG J1021-5835 and GeV J1025- 
5809 and 2EG J1049-5847 and GeV J1047-5840 (see Thomp- 
son et al. 1995, Lamb & MaComb 1997), which are equal to 
~ 10~ 7 ph cm s above 1 GeV. Based on these nor- 
malizations we derive the required acceleration efficiencies 
of nuclei at the shock equal to £ ~ 13.5% (for the monoen- 
ergetic injection) and ~ 12.3% and ~ 8% for the power law 
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log (E / GeV) log (E / GeV) 

Figure 2. (a) The optical depths for 7-rays on pair production in collision with stellar photons (as a function of their energy) 
injected at two distances from the massive star WR 20a and for selected angles measured from the direction defined by the centers of 
the stars, for R = 1.3R WR : a = 49° (thick full curve), 70° (thick dashed), 90° (thick dot-dashed), and 3i? W R: <* = 20° (thin full), 85° 
(thin dashed), 150° (thin dot-dashed), (b) The integral cascade 7-ray spectra escaping from the binary system (marked by the thick 
curves) averaged over the observation angles produced by the primary 7-rays (marked by the thin curves) which in turn are injected 
by protons in collisions with the matter of the stellar winds via decay of pions. Protons are coming from dissociation of nuclei with the 
monoenergetic spectrum (thick full curve) and the power law spectra with the spectral indexes equal to 2. (thick dotted curve) and 2.3 
(thick dashed curve). The cascade 7-ray spectra are normalized to the GeV flux observed from the EGRET sources (GeV J1025-5809 
and GeV J1047-5840, Lamb & MaComb 1997) in the direction of WR 20a The double thin curves show the sensitivity limits for the 
GLAST (1 year observation), MAGIC, and HESS telescopes (5(7 50 hrs observation, > 10 events), e.g. Lorenz 2001. 



injection with the spectral indexes 2 and 2.3, respectively. 
The 7-ray spectra predicted by this model (shown in Fig. 2b) 
should be detectable above 100 GeV by the Cherenkov tele- 
scopes on the Southern hemisphere such as the HESS and 
CANGAROO III or at lower energies by the MAGIC type 
telescope (e.g. HESS II). 

The 7-ray fluxes produced by protons at the shock re- 
gion in collisions with the matter of the winds via decay of 
pions are also accompanied by the high energy neutrinos. Let 
us estimate the muon neutrino event rate from our example 
binary system WR 20a. The spectra of muon neutrinos, cal- 
culated exactly as the primary spectra of 7-rays, are shown 
in Fig. 3a. They are clearly above the atmospheric neutrino 
background (ANB) and above the sensitivity limit of the 1 
km 2 neutrino detector of the IceCube type. In the case of 
monoenergetic injection of nuclei they are also above the 
present sensitivity of the AMANDA II neutrino detector. 
However due to the localization of the source on the South- 
ern hemisphere this model cannot be at present excluded. 
The number of events produced by the muon neutrinos in 
the 1 km 2 neutrino detector of the IceCube type can be 
estimated from 



(9) 



where S = 1 km 2 is the surface of the detector, P V ^ /J ,(E V ) 
is the energy dependent detection probability of muon neu- 
trino (Gaisser & Grillo 1987). 

Applying estimated above values for acceleration effi- 
ciency, the distance to WR 20a equal to 5 kpc, and typical 
wind parameters of the WR stars, M_s = 3, i>3° = 1, we pre- 
dict in the case of the model I ~ 23.5 (23) muon neutrino 
events inside 1 km 2 detector in 1 yr for the case of neutrinos 
arriving from the directions close to the horizon (no absorp- 
tion inside the Earth) and from the nadir (partial absorption 
inside the Earth). The neutrino event rate predicted in the 
case of model II is ~ 21.6 (20.5) and ~ 3.7 (3.5) muon 
neutrino events in 1 km 2 in 1 yr provided that nuclei are ac- 



celerated with the power law spectrum and spectral indexes 
equal to 2 and 2.3, respectively. In these calculations we ap- 
plied the Earth shadowing factors from Gandhi (2000). The 
high event rates predicted for the nuclei injected with the 
flat spectra should be easily tested by the future 0.1 km 2 
Antares telescope. 

A fraction of neutrons, r), dissolved from nuclei can fall 
onto the surfaces of the massive stars. r\ can be calculated 
as a function of the shock distance from the WR star, r, and 
the Lorentz factor of neutrons, 7„, from the formula 



»?(7n) 



2F(cos <9 D ,1) 



^(COS^D, 1) + -F(cOS#D, — COS^d) 



(10) 



where F(x,y) — JJ crA 7 fiWR(e, cos#)(l + cos#)dedcos 6, 60 
is the angle intercepted by the massive star observed from 
the distance r. We have performed such calculations and 
found that r\ weekly depends on the Lorentz factor of nuclei, 
except for the Lorentz factors for which efficient fragmenta- 
tion starts to occur. It is equal to r/ ~ 0.52, 0.33, 0.22 for 
R B h — 1.3i?wR, 2i?wa, and 3i?wR, respectively. The Lorentz 
factors of neutrons extracted from nuclei are equal to the 
Lorentz factors of the parent nuclei 7a = 7n- Neutrons which 
fall onto the surface of the massive stars produce charged pi- 
ons in collisions with the matter of the stellar atmospheres. 
Pions from the first interaction are produced at the charac- 
teristic densities p ~ 3 x 10 14 cm~ 3 (estimated from the p-p 
cross section and applying the characteristic dimension of 
the stellar atmosphere ~ OTJJwr). In such conditions pions 
with the Lorentz factors y n = 10 6 decay before interacting 
with the matter since they pass only ~ 0.3 g cm" 2 . This 
conclusion has been checked by performing the Monte Carlo 
simulations of the interaction process of neutrons with the 
stellar atmosphere which preliminary results has been pub- 
lished in Bartosik, Bednarek & Sierpowska (2003). 

We have calculated the spectra of neutrinos produced 
by neutrons in the atmospheres of the massive stars applying 
the scale break model for hadronic interactions, and taking 
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Figure 3. The differential spectra of muon neutrinos produced 
by protons interacting with the matter of the stellar winds (a) 
and by neutrons interacting with the matter of the stellar atmo- 
sphere (b). Specific neutrino spectra are calculated for the param- 
eters as the corresponding 7-ray spectra shown in Fig. 2b (these 
same curve styles). The dot-dot-dot-dashed curves indicate the 
atmospheric neutrino background (horizontal - upper curve, and 
vertical - lower curve) within a 1° of the source (Lipari 1993), the 
thin dotted line shows the 3 yr sensitivity of the IceCube detec- 
tor (Hill 2001), and the thick dotted line shows the upper limit 
obtained by AMANDA II detector (Ackermann et al. 2004). 



into account the fraction of neutrons falling onto the stars 
77, and the multiple interactions of neutrons with the mat- 
ter of the stellar atmospheres (see Fig. 3b). The neutrino 
event rates, for the acceleration models considered above, 
are estimated on ~ 12.8 (12.6) events for the monoen- 
ergetic spectrum of neutrons, and 11 (10.4) and 2.9 (2.7) 
events for the spectra of neutrons which are close to the 
power law spectra of their parent nuclei with the index 2 
and 2.3, respectively. These event rates are lower than that 
produced by protons since only a fraction of neutrons inter- 
act. From another site this effect is partially compensated 
by the multiple interactions of neutrons with the matter of 
the stellar atmospheres. Note moreover that neutrinos pro- 
duced by neutrons are emitted within the limited range of 
angles, a, around the plane of the binary system defined 
by, tana ~ V-Rmax — Ro/(R — Rwr) (see values defined 
above), i.e within a ~ 79° for 7? max = 2i? WR . This angle is 
much larger than the eclipsing angle of the WR 20a binary 
system ~ 40°. Therefore neutrinos produced by neutrons 
can be also potentially observed from non-eclipsing binaries 
(see Fig. 0. However, in the case of WR 20a, which incli- 
nation angle is ~ 74.5° (Bonanos et al. 2004), the observer 
is located inside the eclipsing cone. The neutrino signal pro- 
duced by neutrons should be modulated with the orbital 
period of the binary system as postulated by the emission 



angle a. This may help extracting the neutrino signal pro- 
duced by neutrons from the atmospheric background. 

The 7-ray and neutrino signals from such close binaries 
as WR 20a can be strongly variable due to the very unsta- 
ble position of the wind collision region (short cooling time 
scales of thermal plasma) . Small change in one of the winds 
may push the wind collision region onto the surface of one of 
the stars changing somewhat conditions in the acceleration 
region. 
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